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Abstract 

While the 12 C(a,7) 16 0 reaction plays a central role in nuclear astrophysics, the cross section at energies relevant to 
hydrostatic helium burning is too small to be directly measured in the laboratory. The /3-delayed a spectrum of 16 N can 
be used to constrain the extrapolation of the El component of the S'-factor; however, with this approach the resulting 
S-factor becomes strongly correlated with the assumed (3a branching ratio. We have remeasured the /3a branching ratio 
by implanting 16 N ions in a segmented Si detector and counting the number of (3a decays relative to the number of 
implantations. Our result, 1.49(5) x 10 -5 , represents a 24% increase compared to the accepted value and implies an 
increase of ss 13 % in the extrapolated S'-factor. 

Keywords: (3 decay, stellar helium burning, properties of specific nuclei 6 < A < 19 
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1. Introduction 

1.1. The astrophysical S-factor of 12 C(a,7) 16 0. 

The rate of the 12 C(a,7) 16 0 reaction, relative to that 
of the triple-a reaction, regulates the relative production 
of carbon (C) and oxygen (O) during hydrostatic helium 
burning in stars, and by doing so has great significance 
for the held of nuclear astrophysics. Not only is the C/O 
ratio at the end of helium burning directly reflected in the 
observed elemental abundances, it also has profound influ¬ 
ence on subsequent stellar evolution: It affects the nucle¬ 
osynthesis of medium-mass and s-process-only nuclei [1], 
long-lived 7-ray emitters [2], and ^-process nuclei [3], it 
affects the chemical composition of white-dwarf stars and 
thus explosion models of type-IA supernovae [4], and it af¬ 
fects the mass of the remnant left behind by core-collapse 
supernovae [5]. 

At the centre-of-mass energies relevant to hydrostatic 
helium burning, E c m . ss 0.3 MeV, the cross section of the 
12 C(a,7) 16 0 reaction is too small to be measured directly 
in the laboratory. Indeed, the lowest point measured so far 
is E c . m . = 0.89 MeV [6]. The prospects of extending these 
measurements down to 0.3 MeV within the foreseeable fu¬ 
ture are incredibly slim, as it would require an improve¬ 
ment of more than 5 orders of magnitude in experimental 
sensitivity. 

A precise and accurate extrapolation to the astrophys- 
ically relevant energies can be made only by including 
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complementary data obtained by indirect methods. At 
0.3 MeV the S-factor of the 12 C(a,7) 16 0 reaction is domi¬ 
nated by El and E2 resonant capture to the ground state 
via the tails of the subthreshold 1“ and 2 + states at 7.12 
and 6.92 MeV, with smaller contributions from cascade 
transitions [7]. 7-ray angular distributions measured at 
higher energies are used to determine the relative contri¬ 
bution of the El and E2 components. 

The extrapolation of the E2 component may be con¬ 
strained using d-wave phase-shift data from (a,a) elastic 
scattering experiments. Using this approach, Tischhauser 
et al. have obtained 8/7.2 (0-3) = keVb [8], which is 

also the value quoted by Buchmann and Barnes in their 
2006 review article [7]. Taking advantage of new (cc, 7) 
capture data, Sayre et al. have recently obtained the more 
precise value of Se 2(0.3) = 62+g keVb [9]. Complemen¬ 
tary constraints on Se 2 (0-3) have been obtained from sub- 
Coulomb a-transfer reactions [10] and a 7-cascade exper¬ 
iment [11]. An alternative approach based on Coulomb 
dissociation has also been explored [12]. 

The extrapolation of the El component is only weakly 
constrained by the p -wave phase-shift data [8, 7[. The 
precision and accuracy of the extrapolation can be signif¬ 
icantly improved by including data from the (3a decay of 
16 N [13]. Using this approach, Azuma et al. have obtained 
S’ei(0.3) = 80(20) keVb [14], which is also the value quoted 
by Buchmann and Barnes [7] . More recently, Tang et al. 
have obtained Sei(0.3) = 84(21) keVb [15]. Other authors 
have used data from sub-Coulomb a-transfer reactions to 
obtain complementary constraints on Sei(0.3) [10]. 

Cascade transitions via the four bound excited states 
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in 16 O were previously believed to contribute anywhere 
between 3 and 32keVb to the total 5(0.3)-factor [7], but 
new data [9, 16] appear to have constrained the their con¬ 
tribution to 11(3) keVb. Combining this with Se 2(0.3) = 
62l® keVb [9] and Sei(0.3) = 84(21) keVb [15], one ob¬ 
tains a total 5-factor of 5(0.3) = 157(23)keVb. Taken 
at face value, the new results of Refs. [9, 16] thus imply a 
significant change in the error budget, with the El ground- 
state capture now making by far the largest contribution 
(13%) to the overall error, while the E2 ground-state cap¬ 
ture and the cascade transitions only make a modest con¬ 
tribution (6%). A precision of 10% has long been de¬ 
sired by astrophysical modellers [17, 18, 19]. This pro¬ 
vides strong motivation for reducing the uncertainty on 
Sei(0.3). 

It should be noted that the uncertainty on the ex¬ 
trapolated 5-factor is a subject of strong debate. Sev¬ 
eral authors have argued that the experimental data is 
compatible with two different values of 5ei(0.3), a high 
value around 80keVb and a low value around lOkeVb, 
see, e.g ., Refs. [20, 21]. Similarly, it has been argued [22] 
that two solutions exist for 5 e 2 (0.3), a high value around 
150keVb and a low value around 60keVb. Furthermore, 
Tang et al. have shown [15] that the 5ei(0.3) value ob¬ 
tained from the simultaneous analysis of phase-shift, cap¬ 
ture and /3ct-decay data is reduced by 20-30% if the old 
phase-shift data of Ref. [23] are replaced with the more 
recent phase-shift data of Ref. [8, 24]. These observa¬ 
tions stand in stark contrast to the very precise value of 
5(0.3) = 161 ± 19(stat)l2( s y s )keVb recently reported by 
Schiirmann et al. [25] based on a global analysis of a se¬ 
lected sample of “world data”. The data selection criteria 
adopted by Schiirmann et al. have since been criticised by 
several authors [22, 26]. 

1.2. The (3-delayed a decay of 16 N 

The focus of the present work is on reducing the uncer¬ 
tainty on iSei( 0.3) by improving the experimental determi¬ 
nation of the /3a decay of 16 N. The only a-decaying state in 
16 0 that is appreciably fed in the (3 decay of 16 N is the 1~ 
state at 9.6 MeV, and as a result this state dominates the 
a spectrum. So much, in fact, that earlier determinations 
of its /3 branching ratio, bp (9.6), have been obtained by 
measuring the total (3a branching ratio, bp a , and assuming 
bp (9.6) = bp a , which is also the approach that we shall fol¬ 
low. There is, however, also a small contribution from the 
high-energy tail of the subthreshold 1“ state at 7.12 MeV, 
and hence a careful measurement of the a spectrum can be 
used to constrain the a width of the 7.12 MeV state, which 
dominates the El component of the ground-state capture 
at the astrophysically relevant energies. The influence of 
the 7.12 MeV state on the a spectrum is enhanced by its 
large (3 branching ratio, bp( 7.12) ~ 5x 10 -2 , relative to the 
(3 branching ratio of the 9.6 MeV state, bp( 9.6) ~ lx 10” 5 . 

While consensus appears to have been established con¬ 
cerning the shape of the experimental a spectrum [15, 27], 


an improved measurement of bp(9.6) is needed to fix the 
absolute normalisation of the spectrum [27]. 

Azuma et al. have shown that the a width of the 
7.12 MeV state, and hence the value of 5ei(0.3), obtained 
from an i?-matrix analysis of the /3-delayed a spectrum 
(simultaneously constrained by elastic-scattering and cap¬ 
ture data) is strongly correlated with the relative feeding 
of the two 1“ states. In particular, a ±9% uncertainty in 
the ratio 6^(9.6)/6^(7.12) results in a ±6keVb uncertainty 
on 5 E i(0.3) [14], 

We present the following, qualitative analysis of how 
a change in the assumed 6^ (9.6) affects the calculated 
Sei(0.3): Consider a level A with a (3 feeding amplitude, 

B\, and a reduced a width, 'yxa- The contribution of that 
level to the total a spectrum is proportional to (Bx^xa) 2 
(following the notation of Barker and War burton [28]). In 
order to scale the spectrum by some factor, /, we must, if 
the general shape of the spectrum is to be preserved, scale 
{BxTXa) 2 for all contributing levels by the same factor /. 

Since Bx for the 7.12 MeV state has been determined ex¬ 
perimentally with high precision [15], any change in the 
contribution of this level to the a spectrum must involve a 
change in its 7 > a . From these considerations we conclude 
that a change in the absolute scale of the a spectrum, i.e. 
bp(9.6) —X fbp( 9.6), must be accompanied by a change 
Txa fTxa f° r the 7.12 MeV state, which, to a reasonable 
approximation, leads to 5ei(0.3) — X /5ei(0.3), consistent 
with the numerical result of Azuma et al. 

We proceed by reviewing what is known experimentally 
about bp( 7.12) and ^(9.6): Tang et al. have recently de¬ 
termined bp( 7.12) = 0.052(2), though details have not yet 
been published. The TUNL evaluation [29] gives bp(9.6) = 
1.20(5) x 10 -5 , implying an uncertainty of 6 % in the ratio 
bp(9.6)/bp(7.12) and hence an uncertainty of ±4keVb in 
5ei (0.3). However, this value of bp (9.6) is based on a single 
measurement that dates back to 1961 [30]. The measure¬ 
ment, performed by the Mainz group, relied on the direct 
counting of /3 particles with a Geiger-Miiller tube. No de¬ 
tails are given on the error estimate, though based on the 
large number of detected a particles (5.0 x 10 4 ) it may 
be concluded that systematic errors dominate. Slightly 
different values can be found in two later publications of 
the Mainz group: Ref. [31] gives ^(9.6) = 1.19(4) x 10 -5 
quoting Ref. [30] as the source, whereas Ref. [32] gives the 
same value, but with an inflated error, bp( 9.6) = 1.19(10) x 10 -5 , 
quoting Ref. [33] as the source, which would appear to be a 
mistake because Ref. [33] makes no mention of the 9.6 MeV 
state, let alone its /3 branching ratio. 

An independent determination of bp (9.6) may be ob¬ 
tained from the relative feeding of the 9.6 and 8.87MeV 
states, determined to be bp(9.6)/bp(8.87) = 1.00(7) x 10~ 3 
by the Mainz group [34], combined with the (3 branching 
ratio of the 8.87 MeV state, given as bp(8.87) = 1.06(7) x 10 -2 
in the current (1993) TUNL evaluation. One thus ob¬ 
tains bp( 9.6) = 1.06(10) x 10 -5 . However, the origin of 
the bp(8.87) value given in the current TUNL evaluation 
is difficult to trace. The 1986 evaluation quotes E. K. War- 
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burton, private communication, as the source, while earlier 
evaluations give a value of 1 . 0 ( 2 ) x 10~ 2 quoting several 
experiments from the 1950s as sources. The reliability of 
the bp( 8.87) value given in the current TUNL evaluation 
is thus difficult for us to assess. 

Recently, Zhao et al. [35] have obtained 6^(9. 6 ) = 1.3(3) 
using an experimental technique similar to the one de¬ 
scribed here. 

1.3. Experimental technique 

We use an experimental technique very different from 
that of the Mainz group to determine bp (9.6). We implant 
a mass-separated, high-energy, 16 N beam in a finely seg¬ 
mented Si detector and measure the two ionisation signals 
produced by the 16 N implantation and the a + 12 C decay 
products that follow with a half-life of 7.13(2) s [29]. This 
allows us to determine bp( 9.6) in a very straightforward 
manner as the number of a decays divided by the num¬ 
ber of implantations, and the difficulties associated with 
absolute counting of /? particles are thus avoided. This 
technique has been used in several previous experiments 
to determine small branching ratios in the /3-delayed par¬ 
ticle decays of 6 He at Lovain-la-Neuve [36, 37], n Li at 
TRIUMF [38], 12 B and 12 N at KVI [39, 40], and 8 B also 
at KVI [41]. In the studies of n Li, 12 B, 12 N and 8 B it was 
possible to correlate the implantation and decay events on 
an individual basis. This was not possible in the study of 
6 He and has not been possible in the present study due to 
the combination of a rather long half-life and a high im¬ 
plantation rate necessary to obtain satisfactory statistical 
precision. 


2. Experiment 

The experiment was carried out at the former Kernfy- 
sisch Versneller Instituut (KVI) in Groningen, The Nether¬ 
lands. A primary 15 N beam was accelerated by the AGOR 
superconducting cyclotron to an energy of 105 MeV and di¬ 
rected onto a CD 2 gas target with a thickness of 6 mg cm~ 2 . 
The secondary beam emerging from the gas target con¬ 
sisted, among other isotopes, of 16 N, produced via {d,p). 
The TRI/rP dual magnetic separator [42] was tuned to se¬ 
lect 16 N 7 + ions with an energy of 80 MeV. At the final 
focal plane 83 % of the secondary beam was identified as 
16 N. 

The detector system, sketched in Fig. 1, consisted of a 
60 pm thick circular Si detector with a diameter of 18 mm, 
and a double-sided Si strip detector (DSSSD) with a thick¬ 
ness of 78 pm and a surface area of 16 mm x 16 mm. The 
two detectors were mounted in a telescope configuration 
with the circular detector serving as A E detector, leaving 
the DSSSD, thick enough to fully stop the 16 N ions, to 
detect their remaining energy. This type of setup provides 
a means to distinguish 16 N from other beam components, 
since the difference in stopping power can be exploited. 
Lastly, two Nal scintillators were placed next to the cham- 
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Fig. 1: Drawing of the detector setup (not to scale). Details on the 
detectors can be found in the text. 


ber to provide 7 -ray identification of 16 N during the initial 
beam tuning. 

48 strips on both sides of the DSSSD, running in per¬ 
pendicular directions, divide the detector into a total of 
2304 pixels, each representing an active volume of approx¬ 
imately 300 pm x 300 pm x 78 pm [43]. The smallness of 
the detection volume implies that the detector is inher¬ 
ently /3 suppressed, i.e., /? particles from the decay of 16 N 
deposit only little energy in one pixel. This has two ad¬ 
vantages: Firstly, the distortion of the a spectrum due to 
/3 summing is minimized, and secondly, the /3-singles spec¬ 
trum does not extend into the energy region relevant for 
the identification of the a-decay branch. 

An a source, consisting of 239 Pu, 241 Am and 244 Cm, 
was used to calibrate the DSSSD. A dynamic range of 0- 
50 MeV is needed for the identification of the implanted 
16 N ions. At the same time, good energy resolution is 
desirable for the measurement of the a spectrum in the 
energy range 1.0-3.5 MeV. To meet both requirements, 
the pre-amplifier signal was split and fed to two amplifier- 
ADC chains with a difference in gain of a factor of 10. The 
low-gain data is used for the identification of the implanted 
16 N ions, while the high-gain data is used for the a decay 
spectroscopy. The energy resolution (FWHM) achieved is 
30keV at 2.4 MeV and 0.9 MeV at 33 MeV. 

The primary beam was operated in on/off mode, with 
the beam gate open for 15 s and then closed for 15 s. A 
logic signal, representing the state of the beam gate, was 
fed to the data acquisition. A clean decay spectrum can 
then be obtained by only including the data collected dur¬ 
ing the beam-off periods. The data acquisition was trig¬ 
gered by a logic OR between signals in the DSSSD, the 
A E detector, and the two Nal detectors. 

3. Data reduction 

3.1. Event reconstruction 

Energy matching of the signals from the front and back 
side of the DSSSD allows efficient suppression of electronic 
noise. It also allows us to disentangle random coincidences, 
which occur with significant probability during beam-on 
periods due to the rather high implantation rate of 10-20 
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kHz. The condition |.Ef ront — ^backl < 2.5 MeV is imposed 
for signals in the low-gain chain and |-Ef ron t — Aback | < 
0.25 MeV for signals in the high-gain chain. 

The DSSSD has a strip width of 300 pm and an in¬ 
terstrip gap of 35 pm, which means that, from purely geo¬ 
metric considerations, 20 % of the detector surface consists 
of interstrip regions (10% on each side). In these regions 
the free charge carriers created by an ionizing particle are 
unlikely to be collected on a single strip, but are instead 
shared between the two strips bordering the interstrip re¬ 
gion. Therefore, when we detect coincident signals from 
adjacent strips, we must consider the possibility that the 
two signals were created by a single particle. If the com¬ 
bined energy matches the energy measured in a strip on 
the opposite side of the detector, we assume that charge 
sharing occurred. 

Our data show that 8% of the a decays suffer from 
charge sharing between front strips and 10 % from charge 
sharing between back strips, in good agreement with the 
geometric estimate. For the 16 N implantations we find that 
5 % of the events suffer from charge sharing between front 
strips while 20 % suffer from charge sharing between back 
strips. A similar front-back asymmetry was also found 
by Torresi et al. in a dedicated study of charge sharing in 
DSSSDs [44], using ions with A-values and energies similar 
to the 16 N ions in the present study. Unlike us, however, 
Torresi et al. find that the probability for sharing between 
front strips is in agreement with the geometric value. Fur¬ 
thermore, they find that a significant fraction of the shar¬ 
ing events on the front side are associated with opposite- 
polarity pulses, implying that they cannot be identified by 
the method used here, which could explain why our value 
is a factor of 2 short of the geometric estimate. We thus 
consider it possible that we fail to identify 5 % of the 16 N 
implantations, and we include this as a systematic uncer¬ 
tainty on our final result. 
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Fig. 2: Energy signal from the A E detector vs. the energy signal 
from the DSSSD. The grey contour shows the cut used to identify 
16 N implantations. 


DSSSD response, but with a AA-signal consistent with 
that expected for 16 N implantations, have also been iden¬ 
tified. The origin of these events, which amount to only 
2 % of the number of identified 16 N implantations, is not 
fully understood. We include them in the final result as a 
systematic uncertainty. 

3.3. Identification of a decays 

With a kinetic energy of 33 MeV the typical implan¬ 
tation depth of the 16 N ions in the DSSSD is 27(2) pm 
with a straggling of 0.33 pm. Since the a particles from 
the decay are fully stopped in less than 10 pm of sili¬ 
con, we can assume full-energy detection for all decays. 
The decay spectrum obtained during beam-off is practi¬ 
cally background-free, see Fig. 3. The high-energy tail of 


3.2. Identification of implanted 16 N 

Using standard methods [45], the energy loss and strag¬ 
gling of 80 MeV 16 N 7 + ions punching through the A E de¬ 
tector are determined to be 47 MeV and 3 MeV, respec¬ 
tively, leaving on average 33 MeV to be deposited in the 
DSSSD. 16 N implantations are identified by applying ap¬ 
propriate cuts to the (F/dsssd? AA)-values, as shown in 
Fig. 2. The vast majority of the 16 N implantations are 
contained in the main locus at -Edsssd ~ 33 MeV and 
A E « 700. Pile-up in the AA-detector produces the ver¬ 
tical band above the main locus. The vertical band below 
the main locus results from 16 N implantations that induce 
only a partial signal in the AA detector because they strike 
the detector close to the periphery of its active area. By 
including these bands we increase the number of identified 
16 N implantations by approximately 5 %. The 16 N implan¬ 
tations constitute 83% of the total number of counts in 
Fig. 2. Other major beam components are 13 C 5+ , 14 N 6+ , 
15 N' + , 13 C 6 + and protons. No other charge state of 16 N is 
implanted in the DSSSD. A class of events with anomalous 
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Fig. 3: Decay spectrum from data taken with the beam gate closed. 

the /3-particle response extends up to 0.8 MeV. Between 
1 MeV and 3.5 MeV we see the signal from the /3a decay of 
16 N. The counts above 3.5 MeV result from the surface of 
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the DSSSD being contaminated with long-lived a-particle 
emitters (originating from standard calibration sources). 
From Fig. 3 we estimate their contribution below 3.5 MeV 
to be negligible. All counts between 1.0 and 3.5 MeV are 
therefore assumed to be 16 N /3a decays. We note that the 
spatial distributions of decay and implantation events are 
in good agreement, as has been verified by a Kolmogorov- 
Smirnov test. 

4. Normalization 

The procedure for counting the 16 N implantations and 
the /3-delayed a decays has been described in Sections 3.2 
and 3.3; however, to determine the branching ratio, bp( 9.6), 
we must also correct for the detection efficiency of the 
experimental setup. Two effects play a main role here: 
Firstly, we only look for decays when the beam gate is 
closed, which means that we must apply a correction fac¬ 
tor, Con/off, that depends on the half-life of 16 N and the 
duration of the beam on/off gates. Secondly, there is the 
issue of the dead time of the data acquisition, which es¬ 
pecially affects the data taken with the beam gate open. 
To determine the necessary correction factors, C d g ad and 
C d ® ad , the raw and accepted trigger signals were moni¬ 
tored. 

4-1. Beam-off 

The total number of decays occurring during the exper¬ 
iment, is related to the number of observed decays, 

^d°e b cV by 

/vrtotal _ /'''foff sy a robs (-\ \ 

decay ^-'dead'-'on/oflr Mecay * 

Defining the duration of the beam-on gate as T and the 
duration of the entire beam cycle as aT, we can write the 
analytical expression for the first correction factor as 

1 — e~ XaT 

Con/off = NT ^gAT _ l )(e~ AT — e~ XaT ) ’ (2) 

where A = log(2)/ti is the decay constant. In deriving 
Eq. (2) we assume the implantation rate to be constant 
over a time period of several beam cycles, which is a good 
approximation in our experiment. Furthermore, we have 
ti = 7.13(2)s [29], T = 15s and a = 2, resulting in the 
correction factor 

Con/off = 2.342(3). (3) 

The second correction factor, C d ® ad , is calculated as the 
ratio between the number of accepted triggers, n D g, and 
the number of raw triggers, N a g, occurring during the 15 
s beam-off period, i.e., C d f ad = N 0 g/n 0 g. The value of 
this correction factor varied throughout the experiment 
between 1.011 and 1.017. 


4-2. Beam-on 

The total number of implantations occurring during 
the experiment, £f, is related to the number of ob¬ 
served implantations, by 

a rtotal rion Arobs (a \ 

iV impl ~ °dead 0 duty Jv impl ' V*) 

The additional correction factor, C duty , has been intro¬ 
duced to account for a 1 kHz chopping of the primary 
beam, necessary for beam intensity monitoring by the cy¬ 
clotron control system. As a consequence of this chopping, 
beam is only delivered during a fraction, D, of the beam- 
on period, determined by the duty cycle of the chopper. 
Let N on be the number of raw triggers and n on the number 
of accepted triggers occurring during the 15 s beam-on pe¬ 
riod, and let N a g and n 0 g denote the same quantities for 
the subsequent 15 s beam-off period. We want to deter¬ 
mine the number of triggers occuring during the fraction, 
D , of the beam-on period where implantations can take 
place. We assume the trigger rates during the remaining 
fraction, (1 — D), of the beam-on period where the beam 
is in fact off, to be equal to the trigger rates during the 
beam-off period. This seems reasonable since the trigger 
rate during the beam-off period is dominated by the 7 - 
ray backgound in the Nal scintillators. Thus, the number 
of raw and accepted triggers occuring during the fraction, 
D, of the beam-on period where implantations actually 
take place, may be obtained by subtracting (1 — D)N a g 
and (1 — D)n 0 g from N on and n on , respectively. We thus 
obtain the following expression, 

Z^on _ Non ~ (1 — D)N Q g 

C dead C duty - non _ {1 _ D)noS - ( 5 ) 

where C d " ad = N on /ri on is the mean dead-time correc¬ 
tion. The convention may seem somewhat artificial, but 
we keep the mean dead-time correction as a separate fac¬ 
tor in order to more easily assess the relative importance 
of the two effects. C d g ad is by far the largest correction, 
varying between 1.26 and 2.61, while the additional cor¬ 
rection, Cduty, varies between 1.025 and 1.087. The large 
variations of C d " ad and, to a lesser degree, Cduty, are due 
to significant variations in the beam intensity. The value 
of D varies between 0.22 and 0.50. To give an idea of the 
relative importance of the various corrections on the final 
result, we list the values averaged over the entire experi¬ 
ment: 

Cf ad = 1.0136(14) 

CT ad = 1.794(7) 

Cduty = 1.054(2), 

where the numbers in parantheses denote the typical sta¬ 
tistical uncertainty for a single 30 s beam cycle, which in 
our view is the most appropriate measure of the overall 
statistical uncertainty on the correction factors. 
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5. Results and discussion 

In total 54 hours of data have been collected, split into 
27 runs. We have identified 1.235 x 10 8 16 N implanta¬ 
tions in the beam-on data and 1467 a decays in the beam- 
off data. The decay spectrum is shown in Fig. 3. The 
maximum is located at an energy of 2.336(6) MeV, where 
no measures have been taken to correct for /3 summing 
or the difference in detector response to a particles and 
12 C ions. Taking into account all the correction factors 
discussed in Section 4, we obtain an a-decay branching 
ratio of bp (9.6) = (1.49 ± 0.04) x 10 -5 . Since the correc¬ 
tion factor, Cdead^duty; is quite large, it is also a poten¬ 
tial source of significant systematic error. In Fig. 4 we 
show the branching ratio obtained for each of the 27 runs. 
The error bars indicate statistical uncertainties. The up¬ 
per panel shows that consistent values are obtained for 
a rather wide range of correction factors (1.3-2.8). The 
lower panel shows that the values are constant with time. 
Thus, we do not find any evidence of systematic errors 
not accounted for in the dead time analysis. The fluctua- 



Fig. 4: Upper panel: The a-decay branching ratio obtained from the 
individual runs shown against the value of defined in 

Eq. (5). Lower panel: The a-decay branching ratio obtained from the 
individual runs. The error bars indicate the statistical uncertainty. 
The solid line shows the weighted mean and the grey area shows the 
uncertainty on this value. 

tions around the average are slightly larger than expected 
from the error bars, resulting in x 2 /dof = 39.3/26. Be¬ 
cause the value of y 2 /dof is not quite satisfactory we find 
it necessary to scale the statistical uncertainty by a factor 
of y / x 2 /dof. Based on the considerations of sections 3.1 
and 3.2 the number of 16 N implantations could be up to 
7 % larger than the number we find, and we include this 
as a systematic uncertainty, giving a final result for the 


branching ratio of 

bp( 9.6)= ^1.49 ± 0.05(stat)/(g ° 0 (sys)j x 10~ 5 . (6) 

This result represents a 24% increase compared to the 
value of 1.20(5) x 10 -5 reported by the Mainz group [30] 
and quoted in the most recent TUNL evaluation [29]. On 
the other hand, it is consistent with the less precise value of 
1.3(3) x 10 -5 recently obtained by Zhao et al. [35] using an 
experimental technique similar to the one described here. 

We have performed an i?-matrix analysis to estimate 
the impact on the extrapolated S-factor [46]. Our analy¬ 
sis suggests an increase of 24% in Sei(0.3) and hence an 
increase of 13 % in the total S'-factor if we adopt the same 
values as in Section 1.1. A similar strong correlation be¬ 
tween bp (9.6) and Sei(0.3) has also been found by other 
authors [14] . It is interesting to note that our new determi¬ 
nation of bp (9.6) results in a slightly improved agreement 
between the preferred Sei(0.3) value implied by the /3a- 
decay data (previously 84keVb, now 104keVb) and the 
Sei (0.3) value implied by the sub-Coulomb a-transfer data 
(lOlkeVb [10]). However, when the error bars are taken 
into account (±21keVb and ±17keVb, respectively) the 
improvement cannot be said to be significant. 

6. Conclusion 

We have measured the (3a branching ratio of 16 N using 
a technique that avoids the difficulties associated with the 
absolute counting of /3 particles. We implant a high-energy 
16 N beam in a finely segmented Si detector and measure 
the energy signals from the implantation and the a decay. 
The branching ratio is then obtained in a straightforward 
manner as the ratio of a decays to implantations. Our 
result, (1.49 ± 0.05(stat)^Q < | ) o (sys)) x 10 -5 , represents a 
24% increase compared to the accepted value [29]. Our 
branching-ratio determination leads to a Sei (0.3) value 
that is 24% higher than previously believed. Adopting the 
same values as in Section 1.1, the total S-factor goes up 
by 13%. Given the significant astrophysical implications 
of such a change in the S-factor, it would be desirable to 
have our result confirmed in an independent experiment. 
Independent confirmation of the /3 branching ratio recently 
reported by Tang et al. for the 7.12MeV state [15], would 
also be of interest. 
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